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INTRODUCTION 





> To introduce a new method of solving real and complex transverse magnetic (TM) and transverse 
electric (TE) pole locations present in Sommerfeld integrand of layered media Green’s functions. 


Purpose: 

> To find out TM and TE poles for single dielectric medium backed by a perfect electric conductor 
topology. 

Scope: 


> Introduce the methodology and explain the differences in formulation for obtaining surface wave 
(SW) and leaky wave (LW) pole locations for both TM and TE characteristic mode equations. 


Plan of development: 


> Develop appropriate polynomial expressions from transcendental equations for SW and LW poles using 
Mittag-Leffler expansion method. 

> Verify the result from ML expansion method against other existing techniques such as por 
approximation technique, fast all modes, and regular perturbation method. 
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> Near field emissions and crosstalk analysis (z = 0). 


> Find TM and TE SW and LW pole locations. 


> Mittag-Leffler Expansion Method. 


> Other techniques involve graphical method, 


Taylor series expansion, regular perturbation 
method. 


> Works for any given: 
> Frequency- f. 
> Real or Complex Permittivity- €, = £p — jE. 
> Electrically thin or thick substrates. 
> Can find more than dominant SW poles. 
> Separate equations for finding SW and LW 
poles 


MISSION 





port? port 4 





port | 


ground plane 


Nanowires on Top/Inside of Microstrip Crossover 
. 2k 
layered media Interconnects’ 
z 
# N+1 
# i+1 sO” < , me 
PE Jez] General representation of stratified 
# i e z" media (for both topologies) a 
# i-1 
g p 
# 0 


“W. Zhou, X. Dai, and C. M. Lieber, “Advances in nanowire bioelectronics,” Reports on Progress in Physics, Vol. 80, 016701 (29pp), 2017. 
7Hoorfar, J. X. Zheng, and D. C. Chang, “Numerical Modeling of Crossover and Other Junction Discontinuities in Two-Layer Microstrip Circuits,” International Journal of Microwave and 


Millimeter-Wave Computer-Aided Engineering, Vol. 2, No. 4, pp. 261-272, April 1992. 


tGolubovic, Ruzica, Athanasios G. Polimeridis, and Juan R. Mosig, "Efficient Algorithms for Computing Sommerfeld Integral Tails," IEEE Transactions on Antennas and 5 


Propagation, vol. 60, no. 5, pp. 2409-2417, Mar 2012. 


Mittag - Leffler Expansion Method 





History and Introduction 


The Mittag-Leffler theorem was introduced by Gösta Mittag-Leffler, a Swedish mathematician. 


Mittag-Leffler expansions are applied to meromorphic functions f(A) whose only singularities are simple poles at A4, A2, À3.... 
With residues (r4, r2, r3.... ) at these poles. In this case, the f(A) can be expressed as follows: 








f(A) = f (0) + ». — + 
Izi n n 


The above expansion is called as Mittag-Leffler expansion and it can be applied to gamma functions, trigonometric 
functions in the complex plane. 


Integro-differential equations containing integrals and/or derivatives, Abel integral equation of second kind and 
general solution of the linear fractional differential equation can all be expressed in terms of Mittag-Leffler 
expansions. 


The major advantage of ML theorem is that it clearly expresses the function f(A) in terms of its poles in complex 
plane, whereas the other conventional methods checks for analyticity and convergence in the region around poles. 


H. J. Haubold, A. M. Mathai, and R. K. Saxena, “Mittag-Leffler Functions and Their Applications,” Journal of Applied Mathematics, vol. 2011, Article ID 298628, 51 pages, 2011. https://doi.org/10.1155/2011/298628. 


Green’s Function for Hertzian Electric Dipole (HED) on a PEC UMKC 
backed Dielectric Medium 
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f €» Hffree space) The Green’s functions contain the TM and TE 
characteristic mode equations, which contains 
P real and complex poles. 


K 
Dry = Ko + j—tan(K,d) 
r 


Drg = Ko — jk, cot(x,d) 


z| Radiation condition is given by: 
X(z=0) 


dl i f $ 3 H (dielectric) A | ki = 2 $ Ski 
Infinite PEC | B 
Fj |$? ki $ > ki 


“S. Barkeshli, “Efficient Approaches for Evaluating the Planar Microstrip Green’s Function and its Applications to the Analysis of Microstrip Antennas,” Ph.D. dissertation, 7 
The Ohio State University, 1988. 





Mittag - Leffler Expansion Method 





e Consider the following substitutions in 


tan(v) = 
Dry and Dr, equations: 


2. (Gm Da E [((2m + = ]4-(2v)? 


v = kd and u = kod 
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2 
5 2) 
vu tan(v) = —jé,-u (TM) oe 
v cot(v) = —ju (TE) e After using first order terms and simplifying them into a 
appropriate polynomial expression for AW poles: 
where 2-2 
ad Epu = (4 FE 1 + : Teos 
u“ =y“ = (SW) in Top sheet of RS E; a (TM) 
cess wate pie 
Tr = ays TF cü 
ut = $- v’ (LW) in Bottom sheet of RS + ( > + 46 re i6 


Invoking the ML expansions for tan(v) and cot(v) — y v+ 4972 ES 26m fa “+ 2a.) 7 *(1412)=0 (TE) 


trigonometric functions [1]: m oyr E ae v= er 


[1]. W. R. LePage, Complex Variables and the Laplace Transform for Engineers. Mineola, NY, USA: Dover Publications, 1961. 





Numerical convergence comparison 


between ML and FAM [2]- tan(v) 
Mittag-Leffler (ML) Expansion 
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Taylor Series (TS) Expansion 





v is real 





8u (—1)"F2°7(2°" — 1) Bon y2n=t 


tan(v) = n)! 


n=1 


No. of terms | Mittag-Leffler | Rel. error % Taylor series | Rel. error % 


expansion (ML) expansion (TS) 

0.05 0.05 1 0.0406 18.925 0.05 0.0834 
0.05 0.05 11 0.0491 1.8394 0.05 0.001 
0.98 1.491 1 1.3006 12.7678 0.98 34.2705 
0.98 1.491 11 1.4729 1.2104 1.4909 0.0027 
1.2 2.5722 1 2.3360 9.1814 1.2 53.3465 
1.2 2.5722 11 2.5501 0.8592 2.5659 0.2430 
1.48 10.9834 1 10.6859 2./086 1.48 86.5251 
1.48 10.9834 11 10.9561 0.2482 8.0998 26.2540 


2. L. Tsang and B. Wu, “Electromagnetic fields of hertzian dipoles in layered media of moderate thickness including the effects of all modes,” IEEE Antennas Wireless Propag. Lett., vol. 6, pp. 316-319, Jul. 2007 


Numerical convergence comparison between ML and TS 
expansions- tan(v) 
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ML Expansion TS Expansion 






v is complex 
Su 


= _—4\n-192n(92Nn _ 
tan(v) = » C an a U)Bon 2n-1 
n=1 


" witogitter | Mer, | Tea Ag 
(ML) 1w 
0.05+j0.02  0.05+j0.02 1 0.0405+j0.0163 18.9308 0.05+j0.02 0.0967 
0.05+j0.02 0.05+j0.02 11 0.0491+j0.0197 1.8397 0.05+j0.02 1.28E-14 
0.98+j0.1  1.444+j0.3143 1 1.2541+j0.2939 12.9422 0.98+j0.1 34.5943 
0.98+j0.1 1.444+j0.3143 11 1.4262+j0.3124 1.2273 1.44434+j0.3143 0.003 
1.2+j0.05  2.5241+j0.3743 1 2.2879+j0.3637 9.2625 1.2+j0.05 53.4238 
1.2+j0.05  2.5241+/j0.3743 11 2.5020+j0.3734 0.8669 2.5210+4j0.3687 0.2476 
1.48+j0.07 6.8776+j5.3490 1 6.5802+/5.3332 3.4177 1.48+j0.07 86.6535 
1.48+j0.07 6.8776+j5.3490 11 6.8503+j5.3477 0.3133 7.2278+j3.0328 26.8859 





ay ad 5 1 
cot(v) = 7 mr ; 
m= = 


1.48 
1.48 
0.98 
0.98 
0.05 
0.02 
0.005 
0.0001 


ML Expansion 


0.091 
0.091 
0.6707 
0.6707 
19.9833 
49.9933 
199.9983 
1E4 


Numerical convergence comparison 


between ML and FAM-cot(v) 


No. of terms 


v is real 


Mittag- 
Kaiia 
expansion 


0.02902 
0.1171 
0.8004 
0.6880 
19.9899 
49.9959 
199.9990 
1E4 


(ML) 


218.7586 
28.6417 
19.3379 


2.5737 
0.0327 
0.0052 
3E-4 
1.3E-7 


sty =) ! 


Rel. error % 


TS Expansion 


Taylor series 


expansion 


0.6757 
0.0910 
1.0204 
0.6707 
20 

50 

200 
1E4 
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(2n)! 


Rel. error % 


(TS) 


642.1202 


0.0012 


52.1386 


2.4E-9 
0.0833 
0.0133 
8.3E-4 
3.3E-7 












1.48+j0.07 
1.48+j0.07 
0.98+j0.1 


0.98+j0.1 


0.05+j0.01 
0.02+j70.01 
0.005+j0.001 


Numerical convergence comparison between ML and 
TS expansions- cot(v) 


ML Expansion 


0.0906-j0.0705 
0.0906-j0.0705 
0.6611-j0.1439 


0.6611-j0.1439 


19.2141-/3.8495 
39.9933-j20.003 
192.31-j38.462 


v is complex 


Mittag-Leffler 


expansion 


0.02897-j0.0605 
0.1167-j0.0692 
0.7907-/0.1303 


0.6784-j0.1421 


19.2206-/3.8482 
39.9959-j20.0020 
192.31-j38.462 


(—1 
cotu = ye n)! 


Rel. error % 


Quia 


173.72 


22.7460 


19.2674 


2.5647 


0.034 
0.0065 
3.3E-4 


TS Expansion 





ONIX 


2 8. 2n- r 


Taylor series 
expansion 


0.6742-j0.0319 
0.0906-j0.0705 
1.0099-j0.1031 


0.6611-j0.1439 

19.2308-j3.8462 
40-j20 

192.31-j38.462 


Rel. error 


% 
1 


509.56 
9E-5 
51.9063 


2.7/E-9 


0.0867 
0.0167 
8.6E-4 


Validation — Pole Approximation Technique [3] 





. f= | GHz, =), d= 1K, Iz-z’ |=0 
Drg = Ko — jk; cot(k;d) —— (a) j : 





Consider, 


e On substituting the p & q transformations in eqn (a) and 
simplifying it further, we get S1 (Real of (a)) & S2 (Imag > 
of (a)). 
e Also consider, p=x+/y. 
x+jy = +1, sin(x + jy) 
e To obtain all types of possible solutions, one should 
consider both signs in above equation. However, to 


demonstrate the method, we shall consider only positive É : ' i f 4 


solution, which is given as: SWP (Approximate Root) 


S1 = x — (L. sin x cosh y — l; cos x sinh y) > Muller’s or Newton-Raphson Method (error tolerance: 107+*).. 








I 1 ; o 2 f(E™) l 
S2 = y — (l; sin x cosh y + L cos x sinh y) E = k2 B (= A ern = Em) — ( d ee) _ g| 
df =) < Tol. 


[3]. M. J. Neve and R. Paknys, “A technique for approximating the location of surface-and leaky-wave poles for a lossy dielectric slab,” IEEE Trans. Antennas and Propagation, vol. 54, no. 1, pp. 115—120, Jan. 
2006. 
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Validation — Pole Approximation Technique [3] UMKC 





a va 
myy 
E & 


i we S1 = x — (L, sin x cosh y — l; cos x sinh y) S2 = y — (l; sin x cosh y + l; cos x sinh y) 
Range of X and Y: -20 < x < 20 and -20 < y < 20 





he N eel meen to AAE 
0 i 
To e A ders final root solutions of (38) 
100 x 100 40.9915 9 0.0048 seconds 
200 x 200 40.7735 5 0.0104 seconds 
500 x 500 40.8277 5 0.0449 seconds 
10.8328 1000 x 1000 40.7831 5 0.1605 seconds 40.8057 


15.2715 seconds 


10000 x 10000 40.8059 4 (+ 20 seconds to 
generate contour graph) 


Second order ML expansion for TM or TE characteristic mode equation takes less than 5 milli seconds to generate the 


polynomial expression and to calculate the initial root. 
[3]. M. J. Neve and R. Paknys, “A technique for approximating the location of surface-and leaky-wave poles for a lossy dielectric slab,” JEEE Trans. Antennas and Propagation, vol. 54, no. 1, pp. 115—120, Jan. 
2006. 
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iam Surface Wave: Initial root to final root convergence using Newton-Raphson UMKC 
eel y method: ML, PAT and FAM methods 





Tolerance = 10714 





No. of 
: l No. of 
No. of iteration : l 
iterations s to final hela 
l, = |kod,/e,—1| | Initial root(ML) Initial root (PAT) Initial root (FAM) | s to final 
to final root a 
root (ML) Ae (FAM) 
TM- 
1.4050 29.2791 5 26.2503 6 32.9678 5 29.8036 
lossless 
TM-lossy 1.4248 28.7321 -j 3.9926 5 27.8984 -j 9.4468 6 32.7026 - j4.9179 5 29.2934 - j3.8288 
TE-lossless 10.8828 40.9742 4 40.8277 5 41.1508 9 40.8057 


TE-lossy 10.9775 41.0119 - j4.2613 4 41.0514 - j4.2800 5 41.3679 - j4.2423 18 41.0226 - j4.25805 


Leaky Wave Convergence to Final Root 


Tolerance = 107 ‘4 


Initial Root 


SU oid elke (Number of iterations to reach final root ) 


Final Root 
parameter 


ML(m=1) ML(m=2) ML(m=3) 


eee 79.4403 — 
e = 
0.9425 >100 27 13 £ j136.1992 


How to identify LWs from ML expansions? 


w- plane location 





0.5242 — j2.1733 


> Leaky waves are infinite in number and it is difficult to find all using polynomial expressions. The graphical 
method is the best available method to identify maximum number of LWs. However, in reality, we need only fixed 


number of LWs that satisfy radiation condition. These can be obtained using ML expansions. 
> Solve for polynomial expression based on the LW substitutions shown in earlier slides. 


> Another way to identify LW poles is by calculating pole location in w-plane 1.e., sine transformation. If = < 


Re(w) < 2 then ıt could be a LW pole. The imaginary of w indicates the behavior of LW pole 1.e., the one that 


satisfies radiation condition or fast-rising ın nature. 
Robert E Collin, "Field theory of guided waves." (1960). 





Numerical versus Asymptotic Evaluation of TE Mode based Sommerfeld Integral 


To find SWP- 
To find LWP- oe 8 -ko + kt e#cot( kz - pa) =0 UMKC 


= | O swp 
f = 1 GHz, e, = 3.25, tan ô = 0, d = 0.14, |z — z'| = 0, l, = |kod Ve, — 1| = 0.9425 (< D1 Iwp located to 
the right of Ky 


Ep = 79.44 — j136.2 Viwp = 0.3669 — [3.2318 
Ciwp = 137.4 + j78.745 

















2 n = 2 X 5 0.4 -y = —- a 
—— Asymptotic os net 
-= = = Numerical 0.2 -=== Numerical || 
1.5 
0 
a — 
Z 1 2 -0.2 
aP, 
i £ 0.4 
0.5 SE 
-0.6 
0 L 
-0.8 
-0.5 — a s — ‘ -1 
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CONCLUSION 


A new method has been successfully tested to obtain TM and TE pole locations present 
in Sommerfeld integrand in the layered media Green’s functions. 





The ML expansion can generate a polynomial form for the transcendental equation by 
incorporating its pole singularities. 


The inverse square term present in tan(v) and cot(v) helps achieve convergence in the 
case of electrically thin substrates. 


In addition, this method can also predict more than one SW pole and several LW poles 
for electrically thick substrates with proper choice of number of expansion terms. 


The future work involves studying convergence issues and limitations of ML expansions 
followed by application to N-layered media Green’s functions. 
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